Spectral weight of doping-induced states in the 2D Hubbard model 
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The spectral weight of states induced in the Mott gap via hole doping in the two-dimensional Hub- 
bard model is studied within cluster dynamical mean field theory combined with finite-temperature 
exact diagonalization. If the cutoff energy is chosen to lie just below the upper Hubbard band, the 
integrated weight per spin is shown to satisfy W+{5) > S {S denotes the total number of holes), in 
agreement with model predictions by Eskes et al. [Phys. Rev. Lett. 67, 1035 (1991)]. However, if 
the cutoff energy is chosen to lie in the range of the pseudogap, W+{5) remains much smaller than 
S and approximately saturates near S « 0.2. ..0.3. The analysis of recent X-ray absorption spec- 
troscopy data therefore depends crucially on the appropriate definition of the integration window. 

PACS. 71. 20. Be Transition metals and alloys - 71.27-|-a Strongly correlated electron systems 



I. INTRODUCTION 

The two-dimensional single-band Hubbard model has 
been widely used to study the role of Coulomb corre- 
lations in the high-Tc cuprates. One of the remarkable 
features of this model was pointed out long ago by Eskes 
et alX, namely, that doping the system with 5 holes does 
not yield unoccupied low-energy states of weight 5/2 per 
spin, like in an ordinary band insulator. Instead, as a 
result of strong local Coulomb interactions, this weight 
is approximately given by W+(5) > 5. The physical rea- 
son for this feature is that both lower and upper Hubbard 
bands must contribute to the generation of itinerant low- 
energy states when holes are added to the system. 

In striking contrast to this prediction, recent X-ray ab- 
sorption spectroscopy (XAS) data on Tl2Ba2Cu06-5 and 
La2-2;Sra;Cu04±5 by Peets et al.^ show a linear behavior 
W+{6) w 5 only up to about S w 0.2. At doping con- 
centrations in the range S — 0.2 . . . 0.3, W+{S) levels off, 
suggesting the inapplicability of the single-band Hubbard 
model for these high-Tc compounds. 

Since the calculations by Eskes et al.^ were carried out 
for small one-dimensional clusters, it is not entirely clear 
to what extent the discrepancies with respect to the data 
in Ref."2. might be related to the simplicity of the theoret- 
ical model. In fact, Phillips and Jarrell^. recently claimed 
that state-of-the-art many-body calculations based on 
the dynamical cluster approximation (DCA)^!^ do indeed 
predict a saturation of W+{6) close to (5 « 0.2, in agree- 
ment with the measurements by Peets et al. 

The aim of this work is to demonstrate that the dop- 
ing variation of the spectral weight of the induced low- 
energy states depends crucially on the choice of the up- 
per limit of the energy window in which these states 
are countedi^ To evaluate W+{5) for the two-dimensional 
Hubbard model we use the cluster extension of dynam- 
ical mean field theory (DMFT)Z.^^ combined with finite- 
temperature exact diagonalization (ED) J!. The results 
show that, if the integration window of W+{6) is cho- 
sen to reach up to the lower edge of the upper Hubbard 
band, then W^+(5) > S, just as predicted by Eskes et alX 
If the cutoff energy, however, is chosen to lie in the range 



of the pseudogap, then W+{5) is much smaller than 5 and 
approximately saturates near 5 « 0.2 . . . 0.3, as found in 
Ref.'^ It is clear, therefore, that the interpretation of the 
XAS measurements must be based on the correct choice 
of the energy window over which the induced low-energy 
states are taken into account. 



II. RESULTS AND DISCUSSION 

The cluster DMFT calculations are carried out for 
the two-dimensional Hubbard model with nearest and 
next-nearest hopping parameters t = 0.25 eV and t' ~ 
-0.075 eV, respectively (band width = 2 eV). The 
onsite Coulomb interaction is U — 2.5 eV and the tem- 
perature is T = 0.01 eV. With this choice of parameters 
the system is a Mott insulator in the zero doping limit. 
To account for intersite correlations, the square lattice is 
viewed as a superlattice consisting of 2 x 2 clusters. De- 
tails of these finite-temperature ED / DMFT calculations 
can be found in Ref<^ 

Figure 1 shows the spectral distributions for a series of 
hole doping concentrations. For simplicity, we show the 
ED cluster spectra as they can be evaluated directly at 
real uj without requiring analytical continuation. Since 
we are primarily concerned here with integrated sections 
of the density of states, these cluster spectra are adequate 
for our analysis. The electron density per spin is given 
by ricr = 0.5(1 — 5). At all dopings, the upper Hubbard 
band is seen to be separated from the low-energy states 
by a broad minimum related to the Mott gap in the half- 
filled limit. The lower edge of this band gradually shifts 
from = 1.5 eV at low doping to about 3.0 eV at large 
doping. The low-energy states close to Ep = reveal 
a pseudogap at about 50 meV for doping up to about 
0.1. In Ref.^ we showed that the origin of this pseudo- 
gap can be traced back to a prominent collective mode 
in the imaginary part of the (tt, 0) component of the self- 
energy. (See also Ref.^) This mode is therefore directly 
linked to spatial fluctuations within the 2x2 clusters and 
cannot be described within single-site DMFT. Electron- 
addition states in the vicinity of this collective mode are 
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FIG. 1: (Color online) Spectral distributions calculated 
within cluster DMFT for two-dimensional Hubbard model at 
several hole dopings (broadening 7 = 0.02 eV, U — 2.5 eV and 
T = 0.01 eV). (a) S = 0.035 . . . 0.20; (b) 5 = 0.31 . . . 0.69. The 
dashed lines denote the bare density of states. The pseudogap 
is located about 50 meV above Ep = at low doping up to 
i5 « 0.1. The lower edge of the upper Hubbard band shifts 
from about 1.5 eV at low doping to 3 eV at large doping. 
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FIG. 2: (Color online) Integrated spectral weight of low- 
energy states per spin as a function of hole doping, (a) 
Integration window extends to minimum below the upper 
Hubbard band, with doping-dependent cutoff energy cjc ~ 
(1 + (5) eV. (b) Integration window extends up to doping- 
independent cutoff energy ojc = 0.1 .. . 3.0 eV. 



highly damped, giving rise to a pseudogap. In fact, the 
energy and strength of this mode exhibit a clear disper- 
sion with doping, which translates into a corresponding 
doping variation of the mean position and width of the 
pseudogap. Moreover, since the collective mode is lo- 
cated slightly above Ep, states in this energy range have 
a much shorter lifetime than states below Ep, giving rise 
to a pronounced particle-hole asymmetry. As shown in 
Refi^ the low-energy region of the DMFT lattice spec- 
tra, obtained via analytical continuation of the cluster 
self-energy to real w, are fully compatible with the clus- 
ter spectra given here in Fig. 1. These lattice spectra 
are also qualitatively consistent with the corresponding 
low-energy range of the DCA spectra for 16-site clusters 
(using t' — 0) derived in Ref.'^ 

The spectra in Fig. 1 demonstrate that the low-energy 
states induced via hole doping are concentrated near Ep 
only at very low doping. With increasing doping, these 
states spread over a larger energy window, until at unit 
doping an uncorrelated empty band appears between Ep 
and the shifted upper band edge at 2 eV. Evidently, upon 
hole doping, spectral weight transfer from the upper Hub- 
bard band is limited to the region close to Ep only at 
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very low doping. At larger doping, unoccupied spectral 
weight is generated in the entire range up to Ep + W. 

Figure 2(a) shows the integrated spectral weight of the 
doping-induced low-energy states, W+{S), where the up- 
per edge uJc of the energy window lies in the broad min- 
imum below the upper Hubbard band. In the present 
case a convenient choice is Wc — {1 + 6) eV. (The low- 
doping behavior of W+((5) was previously shown in Fig. 5 
of Ref j^) The results demonstrate that W+ (6) > 6 in the 
entire doping range, consistent with the predictions by 
Eskes et al} 

Figure 2(b) shows the integrated spectral weight for a 
variety of fixed, doping-independent cutoff energies. For 
UJc in the range of the pseudogap, i.e., Wc ~ 0.2 . . . 0.3 eV, 
W+{S) remains much smaller than 6 and approximately 
saturates near S w 0.2. ..0.3. The reason for this sat- 
uration is that, as pointed out above, spectral weight 
transfer from the upper Hubbard band must proceed at 
larger doping progressively to states farther above Ep, 
eventually covering the range wp to Ep + W . For lUc ~ 
1.0... 1.5 eV, W+{S) > 6, in agreement with Fig. 2(a). 
For LOc > 1.5 eV, the energy window at low doping in- 
cludes part of the upper Hubbard band (see Fig. 1), so 
that W+{S) does not reach zero in the small doping limit. 
In the limit ujc^ (not shown), W+{6) includes all un- 
occupied states. Thus, W+{6) — \ — n„ — 0.5(1 -I- 6). 
To capture properly the integrated spectral weight of the 
doping-induced low-energy states it is evidently neces- 
sary to adjust the cutoff energy so that it roughly tracks 
the lower edge of the upper Hubbard band. 

The cluster DMFT results discussed above are for 
U = 2.5 eV. As can be seen in Fig. 1, the upper Hub- 
bard band is separated from the itinerant unoccupied 
states by a rather wide gap of about 1 ... 2 eV, regard- 
less of doping. The low-energy spectral peak can there- 
fore be defined without ambiguity. Since there is some 
uncertainty concerning the appropriate value of U for the 
high- Tc cuprates, we have calculated the spectral distri- 
butions also for U = 1.5 eV. As shown in Fig. 3(a), the 
upper Hubbard band then lies about 1 eV lower and over- 
laps the upper part of the uncorrelated density of states. 
The distinction between localized and itinerant states is 
therefore less clear than for U = 2.5 eV. Nonntheless, the 
high-energy peak is still separated from the low-energy 
region by a shallow minimum, except in the limit of small 
doping. At this Coulomb energy, the system is barely in- 
sulating at half-filling, so that the separation between 
itinerant and localized states at small doping becomes 
highly ambiguous. 

Figure 3(b) shows the integrated weight of the low- 
energy feature, where the minimum close to the lower 
edge of the upper Hubbard band is used as cutoff en- 
ergy, i.e.. We « 0.7 ... 1.4 eV for 5 = 0.035 . . . 0.32. Be- 
cause of the uncertainties pointed out above, W+{S) does 
not approach zero in the low-doping limit. Neverthe- 
less, W+{d) > S, in agreement with the prediction by 
Eskes et al} Compared to the analogous spectral peak 
for J7 = 2.5 eV, W+{S) is slightly larger for t/ = 1.5 eV. 
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FIG. 3: (Color online) (a) Spectral distributions calculated 
within cluster DMFT for two-dimensional Hubbard model at 
several hole dopings (broadening 7 = 0.02 eV, U = 1.5 eV 
and T = 0.01 eV). The dashed lines denote the bare density 
of states, (b) Integrated spectral weight of low-energy states 
per spin as a function of hole doping. The integration window 
extends to the doping-dependent minimum below the upper 
Hubbard band. Empty dots: U = 1.5 eV; solid dots: U — 
2.5 eV. 



This trend is also consistent with the results discussed in 
Refji The reason for this enhancement is the less clear 
distinction between itinerant and localized unoccupied 
states for smaller U, i.e., some of the incoherent weight 
spreads to lower energies. Inspite of these uncertainties, 
it is evident that W+(<5) does not exhibit saturation. 

The experimental spectra by Peets et al? reveal two 
fairly clearly separated regions: a low-energy feature 
and a peak at about 1.0... 1.5 eV higher energy. It 
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seems plausible therefore to associate the former with the 
doping-induced low-energy states and the latter with the 
upper Hubbard band, in qualitative correspondence with 
the spectra shown in Fig. 1. Since the lower peak has an 
intrinsic width of about 1 eV, it is clearly not possible 
to associate its spectral weight only with the much nar- 
rower pseudogap region. Instead, this low-energy peak 
most likely covers all unoccupied states below the upper 
Hubbard band. The fact that this peak saturates near 
d « 0.2 in the experiment is therefore in conflict with the 
cluster DMFT results for the single-band Hubbard model 
shown in Fig. 2(a). According to the results provided in 
Fig. 2(b), the saturation obtained by Phillips and JarrellS 
appears to be related to the small, doping- independent 
energy. 

III. CONCLUSION 

The integrated spectral weight of states induced in the 
Mott gap of the two-dimensional single-band Hubbard 
model via hole doping is examined within cluster DMFT 



based on finite-temperature exact diagonalization. The 
doping variation of this weight is shown to depend sen- 
sitively on the energy window in which the low-energy 
states are counted. If the cutoff energy is chosen to vary 
with doping and to lie in the minimum below the upper 
Hubbard band, qualitative agreement with model predic- 
tions by Eskes et al. is found. On the other hand, if the 
cutoff energy is taken to be independent of doping, the 
integrated spectral weight exhibits a completely differ- 
ent variation with doping. In particular, at small cutoff 
energies, approximate saturation at low doping is found. 
Since the XAS data by Peets et al. require integration 
over all unoccupied states except for the upper Hubbard 
band, we conclude that these data are not compatible 
with present cluster DMFT predictions for the single- 
band Hubbard model. 
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